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We present a class of very simple inflationary models of two scalar fields which leads to nongaussian 
isothermal perturbations with "blue" spectrum, n > 1. One of the models is inspired by supersym- 
metric theories where light scalar fields naturally acquire masses ~ H during inflation. Another 
model presumes that one of the fields has a nonminimal interaction with gravity £Ra 2 . By a slight 
modification of parameters of these models one can obtain either gaussian isothermal perturbations, 
or nongaussian adiabatic perturbations with n > 1. 
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I. INTRODUCTION 
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It is well known that quantum fluctuations of the in- 
flaton field </> in the simplest inflationary models pro- 
duce gaussian adiabatic perturbations of metric with 
nearly scale independent spectrum, n w 1 OL0]. In the 
models involving several scalar fields one can also ob- 
tain significant isothermal perturbations ||. However, 
such perturbations produce six times greater anisotropy 
of the microwave background radiation than adiabatic 
perturbations (for a review see To obtains small 

MBR anisotropy and still explain galaxy formation from 
isothermal perturbations one should strongly suppress 
entropy perturbations on the horizon scale while keeping 
them sufficiently large on the galaxy scale. This would 
be possible, for example, if entropy perturbations had 
"blue" spectrum decreasing at large scales with n > 1.5. 
However, until now there were no inflationary models 
which would produce such entropy perturbations. In 
principle, it would be possible to obtain such perturba- 
tions in axion models 0] , but in the same theories where 
one could have n > 1.5 one would also have too many 
axion strings with domain walls attached to them 
Therefore typically the authors discussing possible ob- 
servational consequences of inflationary entropy pertur- 
bations concentrated on the models with n < 1, where 
only a minor addition of entropy perturbation to adia- 
batic ones was admissible ||. 

In this paper we will propose a class of very simple 
models where one can get either entropy or adiabatic 
perturbations with n > 1.5. This possibility may be of 
particular interest in relation to open universe models . 
The main idea is almost trivial: If one has a theory of 
a scalar field a with a mass m, then the mode functions 
crfc(i) in the expansion of field operator a(x, i) 
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obey the differential equation 



Here k is a comoving wavenumber. We will normalize 
the scale factor a in such a way that a = 1 at the end of 
inflation. At this time in the most interesting comoving 
wavenumber interval k < H the mode functions |<7fc(t)| 
for to 2 <§; H 2 become Hi 



(3) 



A 2 

which corresponds to the spectral index n « 1 + jjj? > 1. 
The estimate of energy density of fluctuations of the field 
a can be obtained as follows: 
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The main contribution to the energy density is given by 
the modes with ko < k < H, where fco is exponentially 
small, ko ~ H exp(— ^j) ■ Let us consider density per- 
turbations of the field a in the long wavelength limit, 
k <C fco- This limit is most interesting in the models 
with sufficiently large m/H, such that the wavelengths 
fc,^ 1 ~ H~ x exp(^y) after their subsequent growth in an 
expanding universe still remain smaller than the galaxy 
length scale. In this case on the scale fc _1 ^> k$ l the 
fluctuating field a wanders many times in the region 
—H 2 /m < a < H 2 /m, so that its value averaged over 
the domain of a size fc^ 1 vanishes. As a result, addition 
of a perturbation 5a(k) with k <C fco does not lead to 
usual density perturbations m 2 a8(j Perturbations 
of density will be nongaussian and quadratic in 5a; they 
can be very roughly estimated by 
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which implies that 
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In (almost) all realistic models of inflation the ratio 
is extremely small 11 . Therefore initially perturbations 



of the field a practically do not contribute to perturba- 
tions of metric (isocurvature perturbations). However, 
the inflaton field must decay to produce matter. If the 
products of its decay are ultrarelativistic particles, then 
their energy density subsequently decreases as a -4 , where 
a is the scale factor of the universe. Meanwhile, if the 
field a does not decay, or decay very late, then for a long 
time its energy density decreases more slowly, as a~ 3 . As 
a result, perturbations of the field a may give rise to sig- 
nificant perturbations of metric. This is the same mech- 
anism which is responsible for isothermal perturbations 
in axion theory S. However, the axion mass during in- 
flation vanishes. In our case m does not vanish, and eq. 
(^) implies that the perturbations have blue spectrum, 
which is exactly what we need. 

It is not very easy to implement this idea. First of 
all, Hubble constant can be very large in the beginning 
of inflation, and during the last stages of inflation it may 
change significantly. For example, in the chaotic inflation 
scenario in the theory A</> 4 inflation begins at H ~ M p , 
and stops at H ~ vAM p . One must fine-tune the param- 
eters to ensure that H becomes comparable to m exactly 
at the epoch corresponding to the last 60 e- folds, when 
the perturbations responsible for the large scale structure 
of the observable part of the universe have been formed. 
Moreover, prior to this epoch H remains much greater 
than m. But then there is no reason to assume that the 
field a initially was at the point a = 0: it did not have 
any chance to roll down there at H 3> m. 

In this paper we will propose a simple way to over- 
come these problems. It is based on the observation that 
masses of many scalars in supersymmetric theories (ex- 
cept, possibly, the mass of the inflaton field) acquire cor- 
rections Am 2 = aH 2 in the early universe p2[ . Here a is 
some parameter; typically a = O(l), but exceptions are 
possible. To study perturbations of metric which may 
appear in such models we will consider here a simplest 
toy model where the scalar field a acquires the mass ~ H 
during inflation. It describes fields <fr and a with the ef- 
fective potential 
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This potential has two valleys, at a = and at = 0. 
If initially \<f)\ > \a\, the field a rapidly rolls to the valley 
a — and stays there. In this regime the energy density 
during inflation becomes equal to — 



and the Hubble 



constant is H = 
written as follows: 
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Then the potential (^) can be 



V(<j>, a) = ^M 2 cj) 2 + l(m 2 + aH 2 )a 2 
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where a — 4 M i ■ Thus this simple model can be consid- 
ered as a toy model exhibiting generation of an effective 
mass squared aH 2 of the field a during inflation. 

Similar mass term appears during inflation even if the 
fields (j) and a do not interact with each other (g = 0), but 
the field a is nonminimally coupled to gravity. Indeed, 
the term |i?cr 2 does not lead to any corrections to the 
rate of expansion of the universe in the regime a = 0, but 
it gives a contribution 12 fH 2 to the effective mass of the 
field a during inflation fll4|] . In what follows it will not 
be very important which mechanism gives a contribution 
~ H 2 to the mass of the field a; we will just assume that 
during inflation the field a has an effective mass squared 
m 2 (H) = m 2 + aH 2 . Note that if H changes slowly, 
then the field a at large H approaches the state a = as 

follows: a — cr exp (— "'3^* ) ~ o"o ex P ("^ip)- Thus, 
unlike in the situations where the mass is constant, in our 
case the field a within the time ~ (aH) -1 rolls down to 
a = 0, and we may study its quantum fluctuations near 
this point. Note also, that now the relation m 2 (H) = 
aH 2 is satisfied at all large H, so that after fixing the 
parameter a — 0(1) one does not need to fine-tune m to 
be of the same order as H at the end of inflation. Eq. 
1 -Ji suggests that at the end of inflation with aH 2 ^> m 2 
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we will have perturbations Sp ° 
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spectrum, \6%\ ~ (Jr)*" "' Further evolution of these 
perturbations is strongly model-dependent. The main 
point is that under certain conditions the factor 



Ptotal 



may grow up to 1, whereas the shape of the potential may 
not change considerably. In what follows we will study 
one of these possibilities. We will skip the numerical 
coefficients of the order of one since we are interested 
only in order of magnitude estimates. 

According to JTT[ ] , at the stage when the energy den- 
sity is dominated by ^M 2 (j) 2 the scale factor is given 
by a ~ exp(— 2tt4> 2 /M 2 ). Consider the evolution of 
the modes of field a which have the scales bigger than 
the horizon scale H^ 1 during inflation. The comoving 
wavenumber k for these modes satisfies the condition: 
k -c Ha. Neglecting the spatial derivatives in the equa- 
tion for scalar field a and assuming that during infla- 
tion aH 2 m 2 we obtain the following equation for 
(Tfc: 3H&k = —aH 2 (Tk- The solution of this equa- 
tion is <Jk « a~~. At the moment of horizon crossing 
(k ~ Hkdk) the amplitude of quantum fluctuations is 
I <Tfc I k 3 ~ H%, where Hk and afc are correspondingly 
the values of the Hubble constant and scale factor taken 
at the moment when perturbation with the wavenumber 
k crosses the horizon. This together with the relation 
CTfc cx a~~ yields the relation 
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for k <C Ha. Taking into account that Hk 

1 /2 

M (ln(iJfc/fc)) ' , we find that at the end of inflation 
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Here we took the following normalization for the scale 
factor: a = 1 at the end of the inflation when </> ~ 1. 

Detailed evolution of these perturbations after infla- 
tion is model-dependent. For simplicity we assume here 
that the correction aH 2 to the mass of the field a disap- 
pears soon after the end of inflation. This is the case for 
both of our models, but it may not be the case for su- 
persymmetric models studied in fl2|| , where some correc- 
tions to our results will appear. We will also assume that 
the particles a are either stable or decay very late, and 
that the inflaton field decays into ultrarelativistic parti- 
cles with the energy density p r (tf) ~ M 2 immediately 
after inflation. For a long time these particles dominate. 
Hubble constant decays as H oc M/a 2 during radiation 
dominated epoch. The modes with k 3> max {ma, Ha} 
oscillate and their amplitudes decay as a -1 . The modes 
with k -C max {ma, Ha} have practically constant ampli- 
tude at the beginning and then start to oscillate with the 
frequency ~ m when the Hubble constant H(a) drops 
below m. Meanwhile their amplitudes decay as a~ 3 / 2 . 
Using these facts one can easily calculate the spectrum 
of scalar field at a > M/m: 
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To make an estimate of the energy density of fluctu- 
ations of the field a one may calculate the potential" 
energy density 



( P r> = \m 2 (a 2 ) 
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It is easy to check that the contributions of the other 
terms to the total energy either the same order of magni- 
tude or smaller than ([12[). When H drops below m term 
( |l2| ) correctly counts the leading contribution to the en- 
ergy density. Taking into account that <7fc depends only 

and integrating over the angles one can write 



on k 



( [I4 ) in the following manner: 
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Substituting dTTJ) in (|l|) and integrating over k < aM , 
one obtains, for a<Cl, 
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(14) 



The produced u-particles are not distributed homoge- 
neously. If the scalar field is relatively stable and starts 



to dominate in some time after inflation then their en- 
ergy density fluctuations can become very relevant for 
cosmology. As we already mentioned, to estimate these 
fluctuations one cannot simply apply the standard meth- 
ods developed in |l]] because in our case the homogeneous 
component of the field a vanishes. A correlation function 
characterizing energy density fluctuations is 



The power spectrum |<5£| 2 is defined as follows: 
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To make an estimate of |<5£| we will take into account 
only the contribution of the potential term in the energy. 
After simple calculations one finds that 
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Because we are only interested in large scale fluctuations 
which have been produced at inflation, we can make a 
cut off at k ~ M. By substituting the spectrum ( p"l| ) in 
(|l6| ) and taking into account ( |l4| ) one finally gets: 



1^1 ~ C(k) 



(17) 



for k < min(M (f-) 1/2 , Mexp(-l/2a)). Here C(k) = 

^(In^ + ^r^lnf ) 1_f (f Strictly speaking, 
the above expression for \S^\ refers to the relative energy 
density fluctuations in (p a ) after inflation, at a > M/m, 
but before the field a starts to dominate. 

Let us suppose that the particles a at some stage be- 
gan dominating the energy density of the universe. This 
does not necessarily imply that they dominate the en- 
ergy density at present, because they could decay later 
on. However, as soon as they give the dominant contribu- 
tion to density at some moment, perturbations of metric 
at all subsequent moments will be determined by |<S£|. 
Here one should distinguish between two possibilities. If 
the field a does not decay (as in the axion theory) or de- 
cays in a hidden sector without changing the amount of 
photons in the MBR, then the corresponding perturba- 
tions can be called isothermal, or entropy perturbations. 
However, if the field a dominates and decays sufficiently 
early, then we essentially have a secondary reheating, and 
the perturbations (g) should be considered as adiabatic 
perturbations with blue spectrum. In a more general case 
by tuning the moment of decay one can get a mixture of 
adiabatic and isothermal perturbations. 

To evaluate possible implications of this result one 
should first of all estimate the desirable value of the 
coefficient a. Suppose that, as usual, perturbations of 
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metric on the galaxy scale have the wavelength which 
was about 10 25 times greater than H~ x ~ M _1 at the 
end of inflation (our final result will not be very sensi- 
tive to this and the subsequent assumptions). Suppose 
also that C(k) ~ 10 2 for fe ~ 10~ 25 M. Then the re- 
quirement that \5j.\ ~ 10~ 5 on the galaxy scale gives us 
the estimate ^ ~ 0.3, which corresponds to the spec- 
tral index n ~ 1.6. One can easily check that this re- 
sult is consistent with our assumption that C(k) ~ 10 2 
for a broad range of relations between M and m. This 
means that an exact value of a is not very much sensitive 
to our assumptions concerning other parameters of our 
model. (We will specify the possible range of values of 
M and m in Q.) On the other hand, with ?f ~ 0.3 
the amplitude of density perturbations falls more than 
10 times on the way from the galaxy scale to the scale of 
horizon. This strongly suppresses the microwave back- 
ground anisotropy, thus resolving the problem of the 
MBR anisotropy produced by isothermal perturbations 
being 6 times greater than that produced by adiabatic 
perturbations of the same magnitude. Of course it does 
not mean automatically that this model can easily pass 
all cosmological tests, maybe in the end it will be desir- 
able to consider combined effects produced by adiabatic 
and isothermal perturbations. This should be a topic of 
a separate investigation. The main purpose of our work 
is to point out that there exists a very simple but over- 
looked class of models where the perturbations of metric 
have rather unusual properties: They are nongaussian, 
they can be either isothermal or adiabatic, and they have 
blue spectrum with n > 1. 

Note that the nongaussian nature of perturbations in 
these models is especially interesting on scales smaller 



than fe n 



Indeed, on such scales one can introduce an ef- 



fective "classical" quasi-homogeneous field a. This field, 
however, takes different values in different points. The 
calculation of density perturbations in each small region 
of such type can be performed in the standard way |y,^| , 
but the final result will depend on our position in the 
|p|,p^| . This introduces an additional type of the 



universe 



large-scale structure of the universe. For ~ 0.3 the 
corresponding scale ~ fc^ 1 is too small to be of any real 
importance. It would be interesting, however, to study 
the corresponding effects for a<l. 

The models considered above can be easily generalized. 
For example, in supersymmetric theories the terms pro- 
portional to H 2 appear not only in the quadratic parts 
of the effective potential but to other parts as well. A 
more general toy model to study this effect is O] 
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For H <C to ct , the effective mass squared at the minimum 
of the effective potential is 2m 2 . Meanwhile in the e arly 
universe at H 3> m a the minimum is at a — \J ctj (3 M p , 
and the effective mass squared of the field a at that stage 
is aH 2 . If a <C 1, then at the stage of inflation pertur- 



bations of the field a are generated. The main difference 
between this model and the one considered before is that 
in this model during inflation there exists a large ho- 
mogeneous field er. Therefore density perturbations can 
be calculated in a more standard way [lj|^]. They are 
gaussian, with blue spectrum growing at large k as fc"/ 3 . 
In the Affleck-Dine model of baryogenesis [HJ this may 
lead to isothermal perturbations of baryon density with 
the spectral index n ~ 1 + 2a/3. We hope to return to 
this model and other issues discussed above in a separate 
publication p5| . 
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